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Ch , A typical problem with the conventional Galerkin approach for the construction of finite- 

»^ _ mode models is to keep structural properties unaffected in the process of discretization. 

l/^ ' We present two examples of finite-mode approximations that in some respect preserve 

the geometric attributes inherited from their continuous models: a three-component 
model of the barotropic vorticity equation known as Lorenz' maximum simplification 
equations [Tellus, 12, 243-254 (I960)] and a six-component model of the two-dimensional 
Rayleigh-Benard convection problem. It is reviewed that the Lorenz~-1960 model respects 



^ ' both the maximal set of admitted point symmetries and an extension of the noncanonical 

C/3 ' Hamiltonian form (Nambu form). In a similar fashion, it is proved that the famous 

O 



Lorenz~-1963 model violates the structural properties of the Saltzman equations and hence 
cannot be considered as the maximum simplification of the Rayleigh-Benard convection 



C/3 

I problem. Using a six-component truncation, we show that it is again possible retaining both 

symmetries and the Nambu representation in the course of discretization. The conservative 
part of this six-component reduction is related to the Lagrange top equations. Dissipation 
^ . is incorporated using a metric tensor. 

0^ ■ 1 Introduction 

I Various models of the atmospheric sciences are based on nonlinear partial differential equations. 

0^ I Besides numerical simulations of such models, it has been tried over the past fifty years to 

^ ' capture at least some of their characteristic features by deriving reduced and much simplified 

systems of equations. A common way for deriving such reduced models is based on the Galerkin 
approach: One expands the dynamic variables of a model in a truncated Fourier (or some other) 
5^ I series, substitutes this expansion into the governing equations and studies the dynamics of the 

corresponding system of ordinary differential equations for the expansion coefficients. Although 
the number of expansion coefficients is usually minimal to allow for an analytic investigation, 
these reduced models have been used in order to explain some common properties of atmospheric 
models. 

To the best of our knowledge there is up to now no universal criterion for the selection of 
modes or the choice of truncation of the series expansion. However, at least some cornerstones 
for the Galerkin approach are already settled. It is desirable for finite-mode models to retain 
structural properties of the original set of equations, from which they are derived [18, 19]. Such 
properties are, e.g., quadratic nonlinear ities, conservation of energy and one or more vorticity 
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quantities in the nondissipative limit and preservation of the Hamiltonian form. 

Recently, an extension of the Hamiltonian structure based on the idea of Nambu [35] to 
incorporate multiple conserved quantities in a system representation also came into focus. It 
was shown in [7, 36, 38, 40, 43, 47] that various equations of ideal hydrodynamics and magneto- 
hydrodynamics allow for a formal Nambu representation. It therefore seems reasonable to derive 
finite-mode models that also retain this structure. Moreover, almost all models in the atmo- 
spheric sciences possess symmetry properties. These symmetries should thus be taken into 
account in low-dimensional modeling too, which is an issue in the field of equivariant dynamical 
systems (see, e.g., [20]). 

The general motivation for this work is that low-order models are still in widespread use in the 
atmospheric sciences. It has been mentioned above that their original purpose was to identify 
characteristic features of the atmospheric flow in the pre-supercomputer era. While the advent 
of supercomputers partially renders this aim obsolete, finite-mode models are still valuable 
for testing advanced methods in the atmospheric sciences, related to issues of predictability, 
ensemble prediction, data assimilation or stochastic parameterization [2, 32, 41, 42]. Such finite- 
mode models offer the possibility for a conceptual understanding of techniques that are to be 
used in comprehensive atmospheric numerical models later on. For such testing issues, in turn, 
it is essential to have finite-mode models that preserve the structure of the underlying set of 
partial differential equations at least in some minimal way. 

In this paper we give two examples of finite-mode models that retain the above mentioned 
features of their parent model: The first is the three-component Lorenz~1960 model, derived 
as the maximum simplification of the vorticity equation [25]. The second is a six-component 
extension of the Lorenz-1963 model [26]. The authors are aware that there exists a great 
variety of other finite-mode (Lorenz) models, such as e.g. [27, 28, 29], possessing richer geometric 
structure and allowing to address other important issues in the atmospheric sciences, such as the 
existence of a slow manifold, atmospheric attractors, balanced dynamics and the initialization 
problem of numerical weather prediction. Results in these directions can be found, besides in 
the original papers by Lorenz, e.g., in [10, 11, 52, 53]. The choice to investigate the Lorenz- 
1960 and Lorenz-1963 models, however, is reasonable since the latter still is one of the most 
prominent finite-mode models used in dynamic meteorology for testing issues as reviewed above. 
As we are going to show, the Lorenz-1963 model in various respects does not constitute a sound 
geometric model, the derivation of a revised version of this system appears to be well justified. 
The Lorenz-1960 model, on the other hand has been chosen as it is the simplest system for 
which the techniques to be applied in this paper can be demonstrated. 

The Lorenz-1963 model is a dissipative model and as such it necessarily violates conservative 
properties. On the other hand this is a rather typical situation for more comprehensive atmo- 
spheric numerical models too. Usually, the conservative dynamical core of such models is coupled 
to a number of dissipative processes such as friction, precipitation and radiation. Nonetheless, 
it is a necessary condition that the numerics for the dynamical core itself do not violate the 
structural properties of the underlying conservative dynamics [55]. Any valuable toy model of 
the atmosphere should reflect this, e.g. by consisting of the superposition of a conservative part 
and a dissipative part. This is one of the guiding principles for our derivation of the generalized 
Lorenz-1963 model. 
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The organization of the paper is as follows: Properties of discrete and continuous Nambu 
mechanics are briefly reviewed in section 2. Section 3 includes a description of the Lorenz-1960 
model, establishing its Nambu structure and its compatibility with the admitted point symme- 
tries of the barotropic vorticity equation. In section 4, it is shown that the Lorenz-1963 model is 
neither compatible with the corresponding Nambu (Hamilton) form of the Saltzman convection 
equations nor with its point symmetries. We hereafter identify the maximum simplification of 
the Saltzman convection equations [46] that reflects both symmetries and the proper Nambu 
structure of the continuous model. Finally, in section 5 we sum up our results and discuss some 
open questions. 



2 Nambu mechanics 

Since Nambu mechanics emerged from discrete Hamiltonian mechanics, it is convenient to start 
with a short description of the latter. The evolution equation of a general n-dimensional Hamil- 
tonian system is given by 

where F = F{zi) is an arbitrary function of the phase space variables Zj, i = 1, . . . , n, is the 
Hamiltonian function and {.,.} is a Poisson bracket, which satisfies bilinearity, skew-symmetry 
and the Jacobi identity. For discrete Hamiltonian systems, the Poisson bracket is characterized 
by an antisymmetric rank two tensor that can depend on the coordinates of the underlying 
phase space. In modern Hamiltonian dynamics, this tensor is allowed to be singular, leading to 
the notion of a Casimir function C, which Poisson-commutes with all arbitrary functions G{zi) 

{C, G} = 0, V G. 

Setting G = if, it follows that every Casimir is in particular also a conserved quantity. 

Guided by Liouville's theorem stating volume-preservation in phase space, Nambu [35] pro- 
posed a formalism for discrete mechanical systems allowing multiple conserved quantities to 
determine, at the same level of significance, the evolution of a dynamical system. More pre- 
cisely, let us consider a point mechanical system with n degrees of freedom and n — 1 functionally 
independent conserved quantities Hj,j = 1, . . . ,n — 1. The evolution equation for an arbitrary 
function F according to Nambu is 

dF _ d{F,H„H„...,H^_^) _ 

At d{Zl,Z2,...,Zn) 

The above bracket operation is called Nambu bracket, which due to the properties of the Jacobian 
is non-singular, multi- linear and totally antisymmetric. It was demonstrated in [49], that a 
Nambu bracket also fulfills a generalization of the Jacobi identity, which reads 

{{Fi, . . . , Fn+l, . . . , F2n~l} + {Fm {Fi . . . , Fn+2, • • • , F2n-l} 

_|_ . . . _|_ . . . ^ F2n~2, {Fl, . . . , F2„_i}} = {Fi, . . . , {Fn, . . . , i^2n~l}} 

for any set of 2n — 1 functions F^. Various discrete models that allow for a Nambu formulation 
were identified, e.g. the free rigid body [35], a system of three point vortices [36], and the 
conservative Lorenz-1963 model [39], which is discussed in some detail below. 
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It appears that the application of ideas of discrete Nambu mechanics to field equations was 
first considered in [6] (and even earlier in a talk [31]), and later independently by Nevir and 
Blender [38]. It was noted that the singularity of many continuous Poisson brackets of fluid 
mechanics may be formally removed by extending them to tribrackets using explicitly one of 
their Casimir functionals as additional conserved quantity. That is, despite the fact that partial 
differential equations represent systems with infinitely many degrees of freedom, up to now 
there only exist models using one additional conserved quantity. This way, the term continuous 
Nambu mechanics (referring to a Nambu representation of field equations) is at once misleading, 
though it is already used in several papers. 

The restriction to tribrackets may be traced back to the underlying Lie algebras on which 
the Poisson brackets in Eulerian variables are based on [30, 34]. Hence, the fixed relation 
between the dimension of the phase space and the number of conserved quantities used for a 
system representation is lost in continuous Nambu mechanics. In the atmospheric sciences, this 
generalization is called energy-vorticity theory, as the employed Casimir functional is frequently 
related to some vortex integral. Since in the atmospheric sciences the evolution of the rotational 
wind field is dominant over different scales, the energy-vorticity description may be well suited 
for a better understanding of e.g. turbulence. Among others, models that can be cast into 
energy-vorticity form include the inviscid non-divergent 2d and 3d barotropic vorticity equations, 
the quasi-geostrophic potential vorticity equation and the governing equations of ideal fluid 
mechanics as well as equations of magnetohydrodynamics [7, 36, 37, 40, 43]. 

The main problem with continuous Nambu mechanics is that it is up to now not clear whether 
it is possible to state an appropriate condition analog to the generalized Jacobi identity (1) of 
discrete Nambu mechanics. While this is obviously a serious point assessing the self-reliance 
of continuous Nambu mechanics compared to usual noncanonical Hamiltonian field theory, for 
the application of the Nambu formalism this point is not of prior importance. Indeed, the main 
benefit of a continuous Nambu formulation so far lies in the possibility of the construction of con- 
servative numerical schemes. Namely, numerically conserving the antisymmetric Nambu bracket 
automatically leads to a numerical conservation of energy and the second constitutive conserved 
quantity of the bracket. This allows to explain the construction of the celebrated Arakawa 
discretization [3] of the Jacobian operator and enables to generalize the Arakawa method in a 
systematic way to other models possessing a Nambu representation [16, 44, 45, 47]. That is, 
although the Nambu form might appear to be an algebraic curiosity from the theoretical point 
of view, it is nevertheless of high value in the numerical application. This is also the main 
reason, why we aim to care about this form in the course of the present paper. The generalized 
Lorenz-1963 model derived in section 4 is of Nambu form and hence it automatically conserves 
certain representatives of both classes of Casimir functionals of the Rayleigh-Benard convection 
equations. 

It could be argued, that the numerical preservation of only one additional conserved quantity 
is only a little success in view of the infinite number of conserved quantities of two-dimensional 
ideal hydrodynamics. As was demonstrated in [1], not only energy, circulation and enstrophy 
are statistically relevant for the large scale behavior of ideal fluid mechanics, but also is the third 
integrated power of the vorticity. Though this objection can hardly be rebutted, the authors are 
not aware of any truncation conserving an appropriate number of conserved quantities besides 
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the method proposed by Zeithn [54, 55], which however can neither be adopted for ah models in 
fluid mechanics nor in arbitrary geometries. Although similar objections also hold against the 
Nambu bracket approach, the Nambu discretization method nevertheless might be considered as 
an enrichment of existing numerical methods in fluid mechanics (see [14, 48] for a discussion of the 
Nambu discretization in relation to the statistics). Moreover, the Nambu bracket approach goes 
beyond various numerical methods for Hamiltonian field equations, in which solely conservation 
of energy can be assured. 



3 Structural properties of the vorticity equation 

The inviscid barotropic vorticity equation on the /-plane for an incompressible flow can be 
written in form of a conservation law 

| = -[*CI. (2) 

where ip{t, x, y) is the stream function generating two-dimensional nondivergent flow, ( = V'^ip 
is the vorticity and [a,b] = da/dx db/dy — da/dy db/dx denotes the Jacobian. 



3.1 The Nambu structure 

Eqn. (2) possesses an infinite number of conserved quantities, i.e. kinetic energy and all moments 
of vorticity are preserved (see, e.g., [4]). In [38], the Nambu (or energy- vorticity) bracket 



{J-i,J2,^3} :=- / 
Jn 



6Ti 



5C ' 6C 



df, (3) 



was introduced for arbitrary functionals J-i[C],i = 1, . . . ,3. In the above equation, 6/6^ denotes 
the usual variational derivative, df = dxdy is the area element to be integrated within the 
2D-domain ^l. Using appropriate boundary conditions (e.g. cyclic), it can be shown that the 
above bracket is totally antisymmetric. Geometrically, this Nambu bracket is essentially a 
reformulation of the singular Lie-Poisson bracket of ideal fluid mechanics, which is based on 
the infinite-dimensional Lie algebra associated to the group of area preserving diffeomorphisms 

on n. 

Using the bracket (3) it is possible to reformulate eqn. (2) as 

| = (C,£,H}. (4) 

In the above equation, Ti and £ denote the global conserved quantities energy and the second 
moment of vorticity (enstrophy), respectively, which are given by 



^[C] = mfdf, £[C] = \j C'd/. 



3.2 Maximum simplification 

One pioneering work in the field of finite- mode approximations was done by Lorenz [25], who 
introduced a minimal system of hydrodynamic equations based on (2). Using a severe truncation 
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of the Fourier series expansion of (, the fohowing set of ordinary differential equations for the 
remaining three modes A, F, G was derived: 

klAG, (5) 
klAF, 



dF /I 1 



dt \P k^ + P 
dG 1/1 1 



dt 2 VA;2 P 

where k, I are constant wave numbers. We now show, that this model preserves in some sense 
the structure of its continuous counterpart, as the above equations can be derived directly from 
the spectral Nambu bracket of the barotropic vorticity equation. For this purpose, we expand C 
in a double Fourier series on the torus 



C(c) = '^MC' 



Mx 

M 



where M = {mik,m2l) and x = {x,y) are the wavenumber and position vector, respectively. 
Moreover, cm = c^m ~ | (^mi,m2 — iBmi,m2)i where f denotes the complex conjugate. The 
variational derivative may be expanded as 

5C(c) ^^dc^ 

where in the right hand side we consider as a smooth function of cm- Plugging these expres- 
sions into eqn. (3) we find the spectral form of the energy-vorticity bracket 

{Ti,F2,Tz} = y~] kl{mik2 - m2ki)^^—^^ (6) 

^ dCK 5c_(M+K) 5cM 

To simplify the notation, we have assumed that 17 is the unit square. Enstrophy and energy 
in their spectral representations are given by £" = 1/2^j^ckC-k and Ti = 1/2 X^k j^ckc^k, 
respectively. To obtain eqn. (5) from bracket (6), we have to truncate the conserved quantities 
£ and 7i on the set of indices ki,k2 € { — 1, 0, 1} under the following restrictions introduced by 
Lorenz: (i) If the coefficients are real at the onset of evolution they remain real for all times. 
(a) If ci^i = — ci^_i at the onset of evolution this relation holds true for all times. Introducing 
the new variables A = Re(co,i), F = Re(ci^o) and G = Re(ci^_i), it is straightforward to recover 
eqn. (5) from eqn. (6). The maximum simplification of the vorticity equation in Nambu form 
then reads 

^ = kl{V^£ X V^n) =: {z,£,'H} , 
with z = {A, F, G)^, where 



£ = -{A' + F' + 2G'), ?^ = - — + — + 
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2' " 4\P k"^ k^ + l 

The above Nambu bracket is based on the Lie-Poisson bracket of so(3), turning the Lorenz~1960 
model into a particular form of the free rigid body equations. It therefore satisfies all properties 
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of discrete Nambu mechanics. In this respect, the continuous Nambu bracket structure of the 
vorticity equation passes over to the discrete Nambu bracket structure of the free rigid body. 
Note, however, that the Lorenz-1960 model represents a restricted class of a free rigid body, as 
only two moments of inertia are independent. Moreover, it was demonstrated in [8] that the 
above truncation also respects the maximal set of admitted point symmetries in spectral space. 
This set of symmetries consists of discrete mirror transformations {t, x, y, tp) i— t- (t, x, —y, —ip) 
and {t,x,y,ip) i— >■ {t, —x,y, together with combinations of shifts by vr in both x and y 
direction. These shifts are the admitted spectral counterparts of translational symmetries in 
physical space. In particular, using these transformations the above two observations by Lorenz 
can be naturally interpreted as conditions of symmetry. That is, for the Lorenz-1960 model 
preservation of symmetries and preservation of the Nambu structure are mutually compatible. 
Due to the Nambu representation, eqns. (5) also satisfy Liouville's theorem. 

As the Lorenz-1960 model inherits the Nambu structure of the vorticity equation and the se- 
lection of modes can be justified using the admitted point symmetries of the continuous equation, 
the notion of a maximum simplification may be regarded as appropriate. 



4 Structural properties of the Saltzman equations 

In this section, we discuss the structural properties of the convection model derived by Saltz- 
man [46]. That is, we discuss the admitted point symmetries and Nambu form and derive the 
maximum simplification that retains these properties in a minimal form. 

The Saltzman equations we base our investigation on read in nondimensional form [22]: 

As before, ^p is a stream function generating two-dimensional nondivergent flow in the rr-z-plane, 
C = V^V is the vorticity, T is the temperature departure from a linear conduction profile, a is 
the Prandtl number, R is the Rayleigh number and [a, b] := da/dxdb/dz — da/dz db/dx denotes 
the Jacobian. 

In what follows, we aim to distinguish between dissipative and nondissipative systems. In 
the former, we preserve the form of equations as given in (7), while in the latter we neglect 
terms V^ip, V^T. Note, however, that in the second case a different definition of R and a arises, 
see [46] for details. 

Let us consider the domain = [—L,L] x [0,1]. The boundary conditions we adopt are 
free-slip boundaries at both the top and the bottom of the fiuid 

ip{t, x,z = 0) = 7p{t, x,z = l) = 0, C(t, x,z = 0) = C(t, x,z = l) = 0, 

together with 

T{t, x,z = 0)= T{t, x,z = l) = 0. 

Although it could be argued that non-slip boundaries in the vertical would be more natural to 
this viscous problem, the above choice is motivated to be able to incorporate the Lorenz-1963 
model, which is based on free-slip boundaries. In x-direction there are different possibilities, e.g. 
periodic, free-stress or non-slip boundaries. 
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4.1 Symmetries 

We are interested in point symmetries of system (7). For this purpose, let us for the moment 
neglect the impact of the boundary conditions. To compute the maximal Lie invariance algebra 
we used the Maple package DESOLV [12]. The maximal Lie invariance algebra reads 

V = 2tdt + xd, + zd,-i3T-4Rz)dT, dt, 8,, Z{g) = g{t)d^, 

(8a) 

Mf) = f{t)dz + f{t)RdT + f'{t)xd^, X2{h) = h{t)d^ - h'{t)zd^, 

where /, g, h run through the set of real-valued time-dependent functions. Hence, system (7) 
admits scalings, shifts in t and z, respectively, gauging of the stream function and generalized 
Galilean boosts in z- and x-direction, respectively. Moreover, there are two independent discrete 
symmetries given by 

ei: {t,x,z,'ip,T) {t,x,-z,-\lj,-T), 62: {t,x,z,ip,T) 1-^ {t,-x,z,-4!,T). (8b) 

The presence of boundary conditions usually restricts the number of admitted symmetries 
strongly. In the symmetry analysis of differential equations, boundary value problems are rarely 
considered (see [9] for a discussion of this problem). On the other hand, for Rayleigh-Benard 
convection, the consideration of boundaries obviously cannot be omitted. Therefore, we now 
single out those symmetries, which are admitted by the boundary value problem. We are only 
interested in symmetries acting on the space geometry of the problem. This is reasonable as 
transformations acting solely on t, V or T in the course of a series expansion do not place re- 
strictions on the Fourier coefficients and thus cannot be used as a criterion for the selection of 
modes. This at once allows to exclude the transformations generated by dt and 2{g) from our 
considerations. Moreover, scaling generated by V in any case would change the fixed geometry 
in z-direction, so we can exclude it too. The most general transformation generated by the 
remaining basis operators in combination with the discrete symmetries is given by 

{t, X, z, i;, T) ^ {t, 62ix + he2), 6i{z + fei + £3), <5i<52(V' + f'xsi - h'zsi), 6i{T + Rfei)), 

where G M and 6j G { — 1, !}• Acting on the boundaries in z-direction, it is straightforward to 
determine those transformations preserving their values: 

(t, X, z, i/j, T) ^ {t, 62{x + £2), z, 62^^, T), 

(9) 

(t, X, z, i/j, T) ^ {t, 62{x + £2), 1 - 2:, -621P, -T). 

It is now necessary to specify the boundaries in x-direction. A natural choice in the atmospheric 
sciences are periodic boundary conditions. This way, shifts in x-direction are admitted. On the 
other hand, this choice singles out the second discrete symmetry, i.e. we have 62 = 1. This is 
the set of point symmetries on which we subsequently base our truncation. 

4.2 Nambu structure 

To make this paper self-contained, we restate some results given in [7], slightly adapted for the 
special form of (7). The conservative part of system (7) can be represented in continuous Nambu 
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form by 



dC 


'6C 


5n' 




'6C 


sn' 


dt 


6T 


sc. 




TC 




dT 


'6C 


sn 








'dt~~ 


.5C' 











{C,c,n}, 
{T,c,n}, 



(lOa) 



with conserved quantities 



j^{^{V^f -RaTz^df, C = I^C{T-z)df, 



(10b) 



representing the total energy and a circulation- type quantity, respectively. In (10a), the Nambu 
bracket {•,•,•} is defined for arbitrary functional J^j = Ti[C,T],i = 1, . . . ,3 by the equation 



{^1,7^2, ^3} ■- 




+ 



5Fi 



5T ' 6C 



+ 



6C ' 5T 



\df. (10c) 



This Nambu bracket is based on the semi-direct product extension of the Lie algebra of area- 
preserving diffeomorphisms with the vector space of real- valued functions on Q. The process of 
extension of a Lie-Poisson bracket is usually done for systems that incorporate more than one 
field variable. There are different ways how to extend a Lie algebra (see [51] for an excellent 
overview), but the semi-direct extension is common for systems where one variable is advected 
by another (as is the temperature departure by the ^-field). Using the Nambu bracket form, 
the semi-direct product structure can be cast in completely symmetric form. 
We note that there is a second class of Casimir functionals given by 



gdf, 



(lOd) 



where g is an arbitrary function of T — z. This class of Casimirs is not needed for the conservative 
Nambu representation. However, it plays an important role for a geometric incorporation of 
dissipation. Namely, using the generalized free energy 



g = n-s 



^ {V^pf - RaTz - -Ra{T 



df, 



(lOe) 



'nv2^ 2 
where 5 is a realization of the class of Casmirs Sg, it is possible to represent the dissipative 
system (7) via 



^ = {c, c, g} + {c,g), ^ = {T, c, g} + (r, g) , 



(lOf) 



where 



(-^1,-^2) : = 



V-^ + i-^V^^ld/ 



(5C Ra 5T 5T 

is the symmetric bracket of dissipation, which is briefly discussed in the section below. 

A delicate problem is to state appropriate boundary conditions making the Nambu bracket 
(10c) twofold antisymmetric. While the antisymmetry {^i,J^2,-^3} = — {-^ii -^3> -^2} is al- 
ways satisfied due to the properties of the Jacobian operator, the antisymmetry {^1,^^25-^3} = 
— {J-2, J-*!, J-3} follows from an integration by parts. Hence, the boundary conditions must be 
specified in a way such that the resulting boundary terms vanish. For the specified boundary 
conditions in the vertical and periodical boundaries in x-direction, the Nambu bracket is indeed 
completely antisymmetric. 



9 



4.3 Maximum simplification 

In this part, we derive the minimal model of (7) that retains both point symmetries and the 
associated Nambu form of the continuous equations. There are several attempts to extend the 
famous Lorenz-1963 model [26]. The motivation for these extensions is that the Lorenz-1963 
model does not represent characteristic features of Rayleigh-Benard convection properly, as 
noted e.g. in [13]. Several authors have tried to improve the Lorenz-1963 model by attaching 
additional modes, but in various cases this lead to models exhibiting nonphysical behavior such 
as violation of energy or vorticity conservation (e.g. [24]). The problem of energy conservation 
was solved in [50] where a universal criterion for the truncation to energy-conserving finite-mode 
models was established. Moreover, truncation to systems in coupled gyrostats form [18, 19] may 
also lead to models that retain the conservation properties of the original equations. We also 
note, that a single gyrostat is a Nambu system and hence using such a truncation, conservation 
of the underlying geometry may be implemented at least in some minimal form. 

It was shown in [39] that the conservative part of the Lorenz-1963 model allows for a Nambu 
representation via 

dz 

— = Vz-Hi X V,-H2, 
at 

where z = {x,y,z)'^. The conserved quantities are 
ni = -x -az, n2 = -y + -z -rz, 

where r = R/Rc, with Rc being the critical Rayleigh number. However, these two conserved 
quantities are proportional to spectral forms of energy and Jq{T — z)'^df, respectively, while the 
spectral expansion of C under the Lorenz ansatz gives identically zero. Therefore, the Lorenz- 
1963 truncation only allows for a Nambu form that is not directly related to the continuous 
Nambu form presented before. Moreover, if we also try to justify the selection of modes of the 
Lorenz-1963 system using point symmetries, we find that it would be necessary to simultaneously 
use the symmetries 62 and shift in x-direction by 1 , which in any case would violate the boundary 
conditions. That is, the selection of modes is not natural from the symmetry point of view in 
this case. 

Additionally, the Lorenz-1963 truncation does not account for the semi-direct product struc- 
ture of the bracket of the continuous equations. An appropriate discrete realization of this 
semi-direct product structure is given by the special Euclidean algebra se(3) = so(3) x M^. 
The associated Lie-Poisson bracket on the dual se(3)* forms the basis of the Hamiltonian (or 
Nambu) representation of the heavy top equations in the body frame. The Lie-Poisson bracket 
reads [23, 30] 

{F, G} = -U- {VnF X VnC) - T ■ {VnF x VpG + VpF x VnG), 

where 11 and T denote the vectors of angular momentum and the direction of gravity as seen 
from the body, respectively. 

The heavy top model consists of three equations governing the evolution of angular momen- 
tum and three equations for the characterization of the direction of gravity as seen from the 
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body. The maximum simplification of the Saltzman model based on the above Lie-Poisson 
bracket therefore needs a six-component reduction. 

We now proceed with the construction of the modified Lorenz-1963 model. The expansion 
in Fourier series that is compatible with the specified boundary conditions is 



ip = {4>nm sin arnrx + ipnm COS amrx) sin mirz, 

n=Om=l fll) 
oo oo ^ ' 

T = {'&nm sin amrx + 6nm cos amrx) sin mvrz, 



n=0 m=l 

where a is the inverse aspect ratio. 

For the selection of modes for the six-component model, we employ the concept of symmetry 
in a similar fashion as in [8, 13, 22]. In these papers it was demonstrated that the admitted 
point symmetries of the original set of differential equations impose restricting conditions on 
the Fourier expansion that have to be taken into account in the course of the derivation of 
finite-mode models. As discussed above, the equations governing Rayleigh-Benard convection 
admit an infinite-dimensional symmetry group with finite-dimensional subgroup (9) preserving 
the boundary value problem. However, not all symmetry transformations included in (9) may 
be used for a selection of modes since we cannot use all shifts in x-direction as the spectral space 
is essentially discrete. 

For the selection of modes, we aim to use the transformations 

ti: {t,x,z,'il^,T) {t,x + l/a,z,ip,T), 
t2- {t,x,z,'il^,T) {t,x,l - z,-'iIj,-T). 

The task is now to compute the corresponding implications of these transformations on the 
Fourier coefficients, which follow from a straightforward application to the expansion (11). The 
transformation ti e.g. implies 

oo oo 

ip = (0nm sin(an7rx -|- nvr) -|- ipnm cos(an7rx -|- nvr)) sinmvrz, 

n=0 m=l 

and similarly for the transformation of T. Hence ti leads to the spectral transformations 
{(pnm,V>nm,'&nm,dnm) ^ ((-l)"0nm, (-l)"9?nm, (-l)"t?nm, (-l)"6'„m). In a similar fashion, the 
transformation t2 is treated. The corresponding transformations in spectral space hence read 

^1 • {4'nmi ^nmi '&nrm Grim) ' ^ (( 1) 'Pnmj { 1) Vnmi (, 1) ^nmi ( 1) GnmP)i 

m m m „ s i^'^) 

^2 • i'PnmT '^nm i '&nm: Grim) {{~'^) 4'nmi {~^) V^nmi (~1) ^nm; (~1) 



■'nm ) • 



These two transformations give a restriction on the admitted modes since for all n = 2k — 1 
and m = 2/c — 1, respectively, the corresponding Fourier coefficients would violate the symmetry 
property and hence are not allowed in the truncation. 

The discrete symmetry group generated by the transformations (12) is G = {e, ti, ^1^2}) 
where e denotes the identity transformation. By exhaustively studying the implications of 
subgroups of G on truncations of (11), we can derive different low dimensional models in a 
similar fashion as was done in [8] for the barotropic vorticity equation. The list of nontrivial 
subgroups of G is given by Si = {e,ti}, ^2 = {e,t2} and ^3 = {e,t2t-i\. 
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Since we already know that the model to be derived must have six coefficients, it remains to 
select them in accordance with the above subgroups. The first six nonvanishing coefficients under 
consideration of the subgroup are ipoi, Oqi, 4>2i-, '^21, '&21 and ^21- Since these coefficients do 
not include those of the original Lorenz-1963 model, the resulting model will not be considered 
here. Selecting the modes using the subgroup S2 and G itself also does not allow to incorporate 
the Lorenz-1963 model. 

The remaining possibility is given by the subgroup 5*3, leading to the choice of coefficients 
011) V'li) ^^11) ^11) '■P02-, ^02- This choice incorporates both the Lorenz-1963 model and the 
model in [13] and also gives a sound justification for the selection of modes. In addition to the 
symmetries, we also aim to preserve the semi-direct product structure of the Nambu represen- 
tation (10c). For this purpose, the selection of the above listed coefficients based on symmetry 
considerations is still too general. It is necessary to scale these coefficients appropriately. Setting 

<Aii = bA, ipii = bB, V302 = cC, 

^11 = eD, 0u = eE, 602 = fF, 

and plugging the corresponding truncation of (11) into the conservative part of system (7), it is 
found that the scaling coefficients have to satisfy 

_ 1 _ 2a^ 



2b' 7r2(l + a2)' ' 7r262(i + a2)2' 

in order to allow for a Nambu representation of heavy top form: 

{Ti,T2,J'3} ■■= -Vr-Fi • X V^Jg - • {VtF2 X V^Jg + X Vr-Fg). 

Here, tt = {A,B,C)^ and V = {D,E,F)^ are the fluid mechanical analogs of the vector of 
angular momentum and the direction of gravity as seen from the body, respectively. 
The resulting six-component model then reads 

^ "-^ktV^""' 

dC 

Ti7 = 0' 

{eirCE - 2b^fTTBF - 2b^B), 



dt 2be 
dE 

— = -r^fenCD - 2b^f7rAF - 2b^A), 
dt 2be} ^' 

'^ = ^J^{BD-AE). 
dt 2f ^ ' 

The conserved quantities (10b) are correspondingly 

n = -4^((l + a^)b\^{A^ + B^) + 2^3^2 ^ ^Rb'ifaF), 

C = — ^((1 + a^)b^eTiiAD + BE) + 4/7rCF + 4C). 
2ab 
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Using both the heavy top Nambu bracket and the conserved quantities, the finite-mode model (13) 
can be cast in Nambu from 



where x = {A, B,C, D, E, F)'^ . Note, that this Nambu bracket formulation is now structurally 
completely analog to the continuous Nambu bracket (10c) of the convection equations, as it 
is based on the Lie-Poisson bracket using the Lie algebra se(3). In particular, by fixing one 
argument in the heavy top Nambu bracket one naturally recovers all the Hamiltonian properties 
of the Lie-Poisson system, e.g. the Jacobi identity. Bracket (10c) therefore also automatically 
conserves the corresponding truncated forms of the second class of Casimirs (lOd), such as 



It is remarkable that the above model has now a mechanical interpretation similar to the Lorenz- 
1960 model. As the Lorenz-1960 model is a restricted class of the free rigid body equations, 
the finite- mode model (13) in turn may be considered as a restricted class of the heavy top 
equations that is referred to as Lagrange top [23]. This is a heavy top with two moments of 
inertia being equal (a symmetric top) with the position vector of the center-of-mass pointing 
in C-direction (which leads to dC/di = 0). Correspondingly, a number of results valid for the 
Lagrange top might already be passed over to the model (13), such as issues of stability [23] or 
numerical algorithms preserving the Lie-Poisson structure of the above bracket [15]. Moreover, 
due to the additional conserved quantity C = const the Lagrange top is a prominent example 
of a Liouville integrable system [5, 17]. In addition, since dxi/dxi = 0, Vi = 1, . . . , 6, the above 
set of equations also satisfies the Liouville theorem. 

Note that the term linear in C (resp. F) in the expression for C (resp. S) arises due to the 
use of variable T instead of T = T — z and correspondingly do the terms linear in A and B in 
the fourth and fifth equation, respectively. 

If C = at the onset of evolution, the above model reduces to the five-component model given 
in [13] upon rescaling of the Fourier modes. Although in any case C = const during evolution, we 
find it nevertheless important to retain this component in the above model. Firstly, it enables 
to cast the reduced model in Lagrange top form. Secondly, the selection of modes based on 
symmetries does not allow to truncate the Saltzman equations to a five component model, since 
there is no additional criterion that permits one to predict a priori whether to chose (^02 or ^02- 
Hence, both coefficients must be incorporated in the Fourier series expansion. 

The extension to a Nambu-metriplectic finite-mode model is straightforward. Following 
[33], a discrete metric system can defined via dz/dt = gV^P = (z,-P), where z is the phase 
space vector, P denotes a phase space function and g is a tensor. It is further required that 
{Fi,F2) = {F2,Fi) ,\/Fi,F2, which in turn enforces g to be symmetric. In our case, using the 
above truncation and correspondingly the generalized free energy (lOe) as phase space functional, 
we can incorporate the dissipative terms upon using 




S 



Ra 



{3e^Tr{D^ + E'^) + GfirF'^ + UfF) . 



g = -2adiag ( -r^,'2ab'^ 



VT' 



■2(1 + a2) TT\l + a^) 27r2 



) 



Re-^a ' iJeV ' Rp. 



a 
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as metric tensor. Then, attaching y^Q to system (13) gives the maximum simpUfication of the 
dissipative Sahzman equations, which reads 

" ((a^ - ■i)TC-BC + 2eRaE) - (1 + a?)T?aA, 



dt 267r(l + 
dt 



(14) 



^ = ^(evrCi? - 262/vri?F - 2}?B) - (1 + a2)7r2l), 
dt 26e 

^ = -^(evrCZ) - - 26^ A) - (1 + a2)7r2i?, 

dt 2oe 

^ = ^(SZ?-^i^)-4vr2F, 
dt 2/ ^ ^ 

or in the more compact Nambu-metriplectic form 

dx 

— = {x,C,g} + (x,g). 
dt 

System (14) can be considered as a damped Lagrange top and its Nambu-metriplectic form 
completes the geometric picture of the maximal structure-preserving truncation of the dissipative 
Saltzman equations (lOf) discussed in the present paper. 

It was noted in [21] that any vorticity field under the Boussinesq approximation has to satisfy 
the balance equation: 



d_ 

di 



1 2/a _ 2/a _ ^—i 

d 



C dxdz = a 







dz 



(^dx 



z=0 

Straightforward computation for the six-component model shows that this balance equation is 
identically satisfied for any value of the Prandtl number. 

It should be emphasized that the results in this section were derived under the assumption 
of ^2 = 1 in (9). This choice was enforced due to the use of periodic boundary conditions, which 
are of obvious importance in geophysical fluid dynamics. On the other hand, the alternative 
choice of (^2 = —1 would not allow to derive the above discrete convection model, as it would 
require a different selection of Fourier modes (not including the Lorenz-1963 model). How- 
ever, this is quite natural as a physical realization admitting this reflection symmetry involves 
sidewalls at x = —L,L [22]. For such a configuration, the generic Nambu structure (10c) has 
to be supplemented with boundary term contributions, since the second antisymmetry relying 
on integration by parts is then not automatically fulfilled any more. As usual in continuous 
Hamiltonian and Nambu mechanics, the presence of nontrivial boundary conditions complicates 
the appropriate formulation of the models, which in the present case also passes over to the 
finite-mode simplification. 



5 Conclusion and outlook 

In this paper we have addressed the problem of maximum simplification of atmospheric models 
by a discussion of the Lorenz-1960 and Lorenz-1963 model. It was reviewed that the Lorenz- 
1960 model is indeed the maximum simplification of the inviscid barotropic vorticity equation 
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that preserves both point symmetries as well as the inherited Nambu structure of the continuous 
counterpart. This way, we have also implicitly shown that the selection of modes in accordance 
with symmetries is compatible with the inherited Nambu form of the discrete model. Inspired 
from this tutorial example, the Lorenz-1963 model was investigated too. It was found that this 
model neither preserves the proper Nambu structure nor is it compatible with respect to the 
underlying symmetries. This may serve as an additional justification of reported unphysical 
behavior of this model. The proposed extension of the Lorenz-1963 model is a six-component 
truncation that also includes the model [13] as special case. This model is constructed using 
a subgroup of the symmetry group of the Saltzman equations preserving the boundary value 
problem. Moreover, the semi-direct product structure of the Lie-Poisson (or Nambu) bracket of 
the field equations is retained, hence the model is automatically energy- and Casimir-conserving 
in the nondissipative limit. This again implicitly shows that for the presented six-component 
truncation both the Nambu structure and the admitted point symmetries are compatible. In- 
corporation of dissipation leads to a discrete Nambu-metriplectic model, also conserving the 
symmetric structure given by the metric part of the continuous bracket. This compatibility 
of the six-component model with geometric structures may be considered beneficial in view of 
the testing issues reviewed in section 1, for which the newly derived model could be employed. 
Moreover, due to the preservation of important geometric structures, both models that were 
presented in this paper may deserve the notion of a maximum simplification. 

If one aims to use finite-mode models for physical purposes and not only as toy models, the 
question of structure-preserving extensions of such minimal systems of equations is of certain 
interest. It was indicated in section 2 that there is merely the method of Zeitlin that allows to 
construct fully Hamiltonian finite-mode approximations of 2d fluid mechanics. In this method, 
series of approximated n-dimensional models are derived on the two-dimensional torus T^, which 
in the limit n — )• oo converge to the vorticity equation and the equations of a Boussinesq stratified 
fluid, respectively [54, 55]. This convergent sequence of finite-dimensional models exists due to 
the property that for a certain representation of SU{n), in the limit of n — >• oo, the group of 
area-preserving diffeomorphisms (and its semi-direct extension by a vector space, respectively) is 
recovered. One main benefit of this method is that for each finite-mode model a maximal number 
of Casimirs is preserved, which makes such models very attractive, e.g. for the investigation of 
statistical properties of fiuid mechanical systems [1]. For more general settings than or 
for three-dimensional fiuid mechanics, however, the question of a connection between finite- 
and infinite-dimensional Lie-Poisson systems is not fully answered yet. This points to another 
still unsolved question, namely whether it is possible to relate discrete and continuous Nambu 
mechanics in some natural way. Although such a relation was established for two very low 
dimensional models in the present paper, the problem of a proper extension of these minimal 
models without violating the Nambu structure has not been tackled yet. Furthermore, it would 
be interesting to apply the method used in this paper to other models of fiuid mechanics. This 
way, a list of maximal simplified structure-preserving models could be established. This should 
be the issue of forthcoming work. 
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